Thermodynamic scheme of inhomogeneous perfect fluid mixtures 
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Abstract 

We analyze the compatibility between the geometrodynamics and thermodynamics of a binary 

O 

O^' mixture of perfect fluids which describe inhomogeneous cosmological models. We generalize the 

bX). thermodynamic scheme of general relativity to include the chemical potential of the fluid mixture 

. ^ with non- vanishing entropy production. This formalism is then applied to the case of Szekeres 

X' 

;h ' and Stephani families of cosmological models. The compatibility conditions turn out to impose 

symmetry conditions on the cosmological models in such a way that only the limiting case of the 
Friedmann-Robertson- Walker model remains compatible. This result is an additional indication of 
the incompatibility between thermodynamics and relativity. 
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I. INTRODUCTION 



When considering cosmological questions, it is a conventional wisdom to work with the 
homogeneous and isotropic Friedmann-Robertson- Walker (FRW) models, since they repro- 
duce many of the features of the observable universe. Nevertheless, the desire to describe 
details such as the process of structure formation, which leads to local inhomogeneities, 
requires the use of more general models. A natural framework to do this is found in the 
inhomogeneous cosmological models, which are exact solutions to Einstein's field equations, 
representing inhomogeneous universes and which contain the FRW models in a certain geo- 
metric limit, thus, generalizing them. Although the presence of inhomogeneities is crucial in 
the process of structure formation, especially in the context of galaxy formation, the inho- 
mogeneous models to be investigated in this work would be insufficient to take into account 
all the details of the physics of structure formation (see, for instance, \^ for a recent review). 

The most general classes of inhomogeneous models known are the Stephani P| and Szek- 
eres slsolutions, which are irrotational perfect fluid exact solutions and have, in general, no 
isometrics. Both families admit a barotropic equation of state only in the FRW limit, and 
this has risen the question of how to flnd or generate non trivial but physically meaningful 
equations of state. To this question is attached the fact that an evolving universe is, from 
the thermodynamics point of view, a system out of equilibrium, thus, we must make sure 
ihat deflning thermodynamic variables makes sense. In this direction. Coll and Ferrando 
^ have shown that an exact solution admits a "thermodynamic scheme" , that is, has well 
deflned thermodynamic variables compatible with the fleld equations, provided the integra- 
bility conditions of Gibbs equation are satisfled, what turns out to be equivalent to the 
existence of equations of state, not necessarily barotropic. 

Following the thermodynamic scheme proposed in |'4l, Krasihski, Quevedo and Sussman 

nriQ 

[5|, ly, 121] worked out the thermodynamic interpretation of the inhomogeneous cosmological 
models, and found that, in the general case, deflning sound thermodynamics involves forcing 
isometrics on the metrics. However, such analysis was made for a one component perfect 
fluid, complying with particle conservation, which satisfles automatically s = 0, i.e., for a 
flow with null entropy production. 

In an attempt to understand such incompatibility between thermodynamics and perfect 
fluid inhomogeneous cosmological models, in this work we study the simplest non dissipative 
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irreversible case in which the entropy production is in general non vanishing: a binary 
mixture of perfect fluids. To this end, we generalize the thermodynamic scheme to the case 
of a binary mixture, and analyze both the integrability conditions of Gibbs equation and the 
condition of entropy production, due to the fluids mixing process. From the integrability 
conditions we are able to recognize the equations of state which describe the mixture, and 
which do not involve imposing isometrics on the metrics. However, to satisfy the entropy 
production condition the metrics must reduce to those of the FRW cosmologies, a result 
that we interpret as a further indication to the incompatibility between thermodynamics 
and relativity. 

Different aspects of binary mixtures of perfect fluids have been analyzed in general rel- 
ativity. Two-fluid models have been intensively applied to describe cosmological models in 
which one fluid represents radiation and the second one the matter content of the universe 
[8,y,Q,^, 3; 15, 3|- If two- fluids are comoving, then their energy-momentum 

tensor is effectively that of a single fluid. If the fluids are not comoving, then Letelier 
has shown that the energy-momentum tensor can be transformed into a tensor explicitly 
exhibiting a preferred spatial direction which resembles the energy-momentum tensor of a 
viscous fluid. It was also shown in JJ^j and that an anisotropic fluid can be consistently 
described by two-perfect-fluid components. Krisch and Smalley Q] have investigated the 
propagation of discontinuities in the relativistic two-fluid system described by Letelier's ten- 
sor. Zimdahl, Pavon and Maartens |^ have studied inflationary models as a mixture 
of two interacting and reacting fluids within the framework of irreversible thermodynamics. 
The authors consider an out of equilibrium dissipative mixture, with irreversible evolution 
in order to describe the reheating process in inflationary universe models. For this end, t hey 
assume the condition sa = for each component A of the mixture. Finally, Cissoko 22 1 
has analyzed a mixture of a reacting and coupled two perfect fluids flowing with distinct 
4-velocities and introduced the Lorentz factor, associated with the relative velocity of the 
fluids, as an additional thermodynamic variable. In all these works less attention has been 
paid to the study of the compatibility between the thermodynamic and geometric evolution 
of the models, specially when inhomogeneities are present. 

This paper is structured as follows. In Sections II and III we study the thermodynamics 
of a binary mixture of comoving perfect fluids, derive the expression for the entropy pro- 
duction due to the mixing process, and propose a thermodynamic scheme which explicitly 
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includes the chemical potential. In Section IV we apply the conditions for the existence of a 
thermodynamic scheme to the Stephani and Szekeres inhomogeneous cosmological models. 
Finally, in Section V we discuss our results and comment on additional problems for further 
research. 

II. THERMODYNAMICS OF THE SOURCE 

Cosmology addresses, among other questions, the one concerning the process of evolu- 
tion of the universe as a whole. Such process, according to general relativity, involves the 
evolution of both the source and the geometry, connected through Einstein's field equations. 
Besides satisfying the field equations, the source must be subject to the laws of thermody- 
namics, which formally, involve quantities defined only in equilibrium. Thus, to describe 
an evolving system, from the thermodynamic viewpoint, we must extend the concepts and 
laws to the non equilibrium case. The simplest way to do this, and the most invoked one 
in cosmology (implicitly or explicitly), is based on the local equilibrium hypothesis, which 
asserts that every point of a system not too far from equilibrium, has a neighborhood which 
is in equilibrium, that is, in such a neighborhood we can define uniquely all thermodynamic 
variables, with the same physical meanings as in equilibrium and satisfying all the relations 
the variables satisfy in equilibrium. Thus, the gradients and temporal changes of the vari- 
ables in the system are connected through the relations they satisfy in equilibrium, i.e., we 
must treat them as scalar fields connected through the equilibrium relations. With this in 
mind, we now consider the source of the field in our models: a perfect fiuid consisting of a 
binary mixture of perfect fiuids. Its energy-momentum tensor is 

TaH = (P + P)UaUl3 + PQafS (l) 

where p and p are the total matter-energy density and the mixture's pressure, respectively, 
and Ua, (« = 0, 1, 2, 3) is the 4- velocity associated with the total matter fiow. Such energy- 
momentum tensor corresponds to a continuum in which energy transport occurs only due 
to matter transport. Thus, entropy production will be a consequence of the mixing process 
alone, as the composition of the mixture will be, in general, inhomogeneous. This tensor 
satisfies the conservation law T'^^ — 0, which in particular, implies the contracted Bianchi 
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identity 

p+{p + p)Q = (2) 

where B = u°'.^ is the scalar expansion. Besides, we restrict ourselves to the case in which 
particles of both components are conserved 

= (3) 

where n is the total particle number density. The analogues of equations Q and (jHl) for 
the case of a one component perfect fluid, together with the corresponding Gibbs equation 
imply, automatically, null entropy production, s = 0. However, for the case of a binary 
mixture, Gibbs equation reads 

Tds = d{p/n) + pd{l/n) — ^idc — fi2dc2 (4) 

where s is the entropy per particle, T is the temperature, pi and /i2 are the chemical 
potentials of components 1 and 2, and c and C2 are the fractional concentrations (c + C2 = 1) 
of each kind of particle in the fluid element under consideration. It is clear that the fluid 
composition is determined by only one parameter, so we can write 

Tds = d{p/n) + pd{l/n) — fide (5) 

where we have introduced the mixture's "chemical potential", p = pi — p2- Equation Q, 
together with Q and (jHl), now implies 

s = -^c (6) 

This corresponds to the entropy production in a fluid in which the mixing process is the only 
irreversible process taking place, and it depends only on the change of the Gibbs free energy, 
due to change in composition, and on temperature. This represents the main difference with 
respect to the case of a one component fluid, and it is our interest to exhibit its consequences. 

III. GIBBS 1-FORM 

As equation ® states, to describe the thermodynamics of any fluid element of the mix- 
ture, we must introduce the conjugate pair of variables p and c, in terms of which the entropy 
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production is written. Thus, these variables play an important role in the description of ther- 
modynamic evolution. On the other hand, the field equations determine p and p and their 
evolution (through the conservation law T^^^.^ = 0), linking geometrical and physical evolu- 
tion. What conditions are necessary and sufficient to guarantee compatibility between the 
field equations and the laws of thermodynamics? According to Coll and Ferrando 4], it is 
enough to demand the compliance of the integrability conditions of Gibbs equation. Thus 
we consider the Gibbs 1-form 

n = Ts „dx" = [(p/n),„ + - /ic „]dx" (7) 

where a comma denotes partial differentiation and d is the exterior derivative. Its integra- 
bility conditions are 

dQ = (8) 

dfi A = (9) 

the ffrst of which is a sufficient condition, of no physical significance, as it involves the exis- 
tence of a "heat function" which would imply unphysical properties to the fiuid. The second 
condition, however, is a necessary and sufficient condition. It represents the differential re- 
lations that the equations of state (4 in the case of a binary mixture in 1 phase) must satisfy 
in order to contain all the equilibrium information about the system, since not all of them 
are independent (only 3 in our case), as one of them can be deduced from the others via the 



Gibbs-Duhem equation. Its components are 

dnAQ = Ztijdt A dx' A dx^ + %fedx* A dx^ A dx^ (10) 

Ztij = -T (^-^pitn^i + /i[,tc sj] = (11) 

Zijk = -T (^Pl^n,j + PiiCj) = (12) 



where square brackets denote antisymmetrization. Using Gibbs 1-form, one can easily see 
that for the case c = 1, i.e., for a one component ffuid, one obtains conditions equivalent to 
those used by Krasihski et al. 0,0, Q]- 

Then, we say that a given cosmological model with a binary mixture of perfect ffuids as 
source satisfies the thermodynamic scheme if the integrability conditions (lllj) and ()12|) are 
fulfilled and if it evolves in accordance with the entropy production law ©. 
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IV. EQUATIONS OF STATE AND EVOLUTION 



In this section we exploit the necessary and sufficient integrabihty condition, derived in 
the last section, in order to find the relations linking the new variables fi and c with the old 
ones, and establish the consequences of the entropy production condition. We do this for 
two general families of inhomogeneous solutions of Einstein's equations, namely the Szekeres 
and the Stephani solutions. 

A. Szekeres solutions 

The source in these models is an irrotational and geodesic perfect fiuid, such that we can 
find local comoving coordinates in terms of which the metric has the form 2^ 

els'" = de - e^'^dz^ - e^%dx^ + dy^) (13) 

with a = a(t,x,y, z), P = Pit, x, y, z), functions to be determined from the field equations. 
In these coordinates the 4- velocity is m'^ = 6q, which implies ii'^ = and so p = p{t). For 
their study, we must consider separately the cases /3' = and /3 7^ 0, where P' = dp/dz. 
In the case P' = 0, the solution is given by 

= 5 (14) 

= A + <I>S (15) 

where $ = $(t), k is an arbitrary constant, A = X{t, z), S is determined by 

^ _ |?7(x^ + y') + V^x + V^y + 2W 

l + lk{x' + y2) ^'^^ 

with U = U{z), V\ = Vi{z), V2 = V2{z), W = W{z) arbitrary functions, $ is given by 
and A satisfies 

X^tt^ + A,t$,t + X^^tt + X^np = U + kW (18) 
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Both equations can be solved once the choice p = p{t) has been made. The matter-energy 
density is given by 



This family of spacetimes has in general no isometries, but when (A/$)^j = the solution 
reduces to a FRW model. 
In the case 7^ we have 



= (20) 



he-'^ie"),, (21) 

with $ = $(t, z), u = z/(x, y, z), h = h{z) and 

e-" = A{x^ + y^) + 2Bxx + 2B2y + C (22) 

where A = v4(z), i?i = -Bi(z), -B2 = B2{z), C = C{z) and h{z) are arbitrary functions, $ is 
defined by 

where p = p{t) is an arbitrary function and k = k{z) must satisfy 

AC - Bl -Bl='^ {h-^ + k) (24) 

This family has in general no isometries, but the FRW model results when $ = zR{t) and 
k = koz"^, where fco is a constant. We now analyze for both families the conditions under 
which the thermodynamic scheme is satisfied. 

Let us consider condition Zij^ = for this solution, in which p,i = 0, and so according to 
Eq.dT^ it reduces to 

fJ'liCjS^k] = (25) 

If we interpret this condition as an algebraic equation for the spatial gradients, V, of the 
variables entering it, the general solution of this equation can be written as Vs = aV/x+feVc, 
with a and b arbitrary functions depending, in general, on all coordinates (t,x*). On the 
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other hand, the spatial components of the Gibbs 1-form (|7j) for a geodesic fluid imply the 
relation 

Hence we can identify a = 1/T and b = —fi/T so that V/i = V[(p + p)/n]. The integration 
of this last equation 

fi = 27) 

n 

yields our first equation of state where we are neglecting an additive function of time for 
two reasons: first, in a two-component system in one phase we have only three degrees of 
freedom, which, once fixed, determine all the remaining variables in such a way that if we 
consider p, p and n as independent, we cannot add to the variable p an arbitrary function of 
time. Second, a chemical potential of this form can be interpreted in a natural way. In fact, 
considering the Gibbs 1-form as given in Eq.Q and the chemical potential (j27jl . we have 

ds = 4^dp - -^dn - ^dc (28) 
riT ^ nT T ^ ^ 

i.e., the chemical potential of the fluid mixture must determine both the matter flow due 

to inhomogeneities in the total particle number density (with uniform composition) and the 

matter flow due to inhomogeneities in the composition (with uniform total particle number 

density). Now, if we substitute p in the spatial components of the Gibbs 1-form, in terms 

of p and p, which are determined from the field equations, and of n (remember that particle 

conservation implies n = f/y/A, where /(x*) is an arbitrary function of spatial coordinates 

and A = det{gij)), then we are in position to recognize the rest of the equations of state of 

the mixture 

c = ln^ T=^ . = ^ (29) 

n n T 

where Cq and Tq are functions of time only. Thus, the system can accommodate, among 
others, the ideal gas equation of state, which, as simple as it is, is more realistic than the 
barotropic one. 

Consider now the condition Zuj = 0. For a geodesic fluid it can be written as the equation 
-?r Vs xVn + sVc xVu + iiVs x Vc + cVu x Vs = (30) 
which by using the components of the Gibbs 1-form ((Tj) reduces to 

- Vs X f- Vn + Vc) = (31) 



n \n 
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This condition is trivially satisfied in view of the equation of state for c. Finally, we must 
demand that the set of thermodynamic variables satisfy the entropy production condition 
(jH)). By inserting the corresponding values given in Eqs.(|27|) and (j29|) . it can be shown that all 
the thermodynamic variables become functions of time only, and that the arbitrary functions 
which define the metric take the values corresponding to the FRW limit. Consequently, the 
Szekeres family of solutions does not satisfy the thermodynamic scheme in general, but only 
in the limiting FRW case. 



B. Stephani solutions 



These models represent the most general conformally flat solution with an irrotational 
perfect fluid as source, admit in general no isometrics and generalize the FRW models. We 



can find local comoving coordinates in terms of which the metric has the form 



23| 



where 



D 



FVt 
V 



h 

4 



(32) 

(33) 
(34) 



(x - xo)^ + {y- yof + {.z- zqY 

and F{t), R{t), k{t), Xo{t), yo{t) and ZQ{t) are arbitrary functions of time. The 4-velocity 
and state variables for this metric are given by 



u 



p = -3C2 + 2 



VCC 



n = fV^ 

where /(x*) is an arbitrary function of spatial coordinates and C{t) is defined by 



(35) 
(36) 
(37) 
(38) 



k = R' 



F2 



(39) 



The FRW limit for this solution is obtained when k, xq, yo and zq are all constants, or 
equivalently, when Vt/V is independent of the spatial coordinates {x, y, z}. 
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To study the thermodynamic behavior of these solutions as a binary mixture of perfect 
fluids, we first analyze the condition Zijk = 0. Then, considering that in this case the 
thermodynamic system is the same as in the Szekeres models, we demand that the equations 
of state [cf. Eqs.(jZ3) and ^] 

H = s = - (40) 

n I 

be satisfied. Moreover, introducing Eqs. ()40|) into the Gibbs 1-form and considering that now 
p^i = 0, we find that 

c = In ^ (41) 

ns 

with Co again a function of time only. On the other hand, the condition Zuj = can be 
written as 

4 Vp X (-Vn + Vc) =0 (42) 

where we have used Eq. ()40|) and the components of the Gibbs 1-form. Considering now the 
explicit form of the equations of state for the Stephani solutions, it can be shown that the 
integrability condition (jl^ is satisfied if 

r^i^H. (43) 

where To is a function of time only. Furthermore, we have to demand that the thermody- 
namic variables so defined satisfy the entropy production condition (jHl). This can be shown 
to be equivalent to demanding that the expression Vt/V depends on time only, a condition 
that reduces the Stephani metric to its FRW limit. 



V. CONCLUDING REMARKS 



We have introduced the set of thermodynamic variables to describe the most general 
inhomogeneous cosmological models, considering that the source is a binary mixture of per- 
fect fluids with, in general, non vanishing entropy production. This process does not involve 
forcing isometrics on the metrics. Nevertheless, demanding the compliance of the entropy 
production condition reduces the metrics to those of the FRW models. This condition must 
be regarded as a consequence of the local equilibrium hypothesis, so our results indicate 
certain incompatibility between thermodynamics and relativity. 
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We believe that in order to generalize thermodynamics to the relativistic case, one should 
first answer questions about how should one incorporate the second law of thermodynamics 
into the thermodynamic scheme and about the invariance of the thermodynamic variables. 
In particular, we note that when irreversible processes appear due to inhomogeneities in a 
system, their main effect is, precisely, the vanishing of the inhomogeneities, indicating the 
need to have a different thermodynamic framework, in which non uniform equilibrium states 
should be possible. 
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